Dynamics of two interacting electrons in a one-dimensional crystal with impurities 
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We investigated the role that the electron-electron interaction plays on the propagating properties 
of wave packets in a one-dimensional crystal with impurities. We considered two interacting particles 
with opposite spins in a band, where we treated their interaction along the Hubbard model. We 
have obtained the density of states of the crystal for different values of the interaction term, as well 
as solved the dynamical Schrodinger equation by varying the initial conditions. We have introduced 
a method through which we were able to follow the time evolution of the wave packets for both spins 
showed in three-dimensional plots, and have evaluated, for each particle, the corresponding MSD's 
and the centroids as function of time . These measurements allow us to determine the influence of 
the interaction on dynamical properties. We discussed the combined effect that the extension of 
the initial wave packets and the interaction strength have on propagating properties. Under certain 
conditions we obtained an entanglement of the two packets associated with both spins that takes 
place in a small region of the lattice. 

PACS numbers: 72.15.Rn, 73.20.Jc, 78.30.Ly 



I. INTRODUCTION 

According to the scaling theory of localization devel- 
oped by Abrahams et ali in low dimensional systems (one 
and two dimensions), any degree of disorder will prevent 
the appearance of a metallic phase. Moreover, former 
experiments done with two-dimensional electron systems 
fabricated on semiconductor surfaces showed a logarith- 
mic increase of resistivity while lowering the tempera- 
tur o^i'^i'^1^ . This behavior lends support to the scaling 
theory of localization, since this happens in case of a weak 
electron-electron interaction. These results were in agree- 
ment with theoretical predictions^' that weak electron- 
electron interaction increases localization. The above 
experiments were done with samples of high density of 
electrons, i.e. systems for which the associated Wigner 
radius Ts 1, which is the ratio between the Coulomb 
energy to Fermi energy. 

However, such a scaling theory does not take into ac- 
count the electron-electron interaction which was lately 
believed to be responsible for a metal-insulator transi- 
tion, observed in several experiments performed at zero 
magnetic field described below. In the metallic phase, 
one observes a strong temperature dependence (a steep 
dg/dT > 0) caused by the delocahzing effects produced 
by the interaction between the particles. 

As the density is reduced such that ^ 1, the inter- 
action becomes dominant, for that regime Finkelstein^ 
and Castellani et ali^ predicted that for sufficient strong 
interactions, a 2D system should present a conducting 
phase as the temperature is lowered. Since recently was 
possible the fabrication of 2D samples of high quality 
with very small amount of randomness, measurements 
were done at very low particles densities. In this way the 
strongly interacting regime ( Ts 3> 1) has become experi- 
mentally accessible. For instance for Ts > 10 experiments 



done on low-disordered 2D silicon samples demonstrated 
that with increasing electron density one can cross from 
the insulating regime where the resistance diverges with 
decreasing temperature, to a regime where the resistance 
decreases strongly with decreasing temperature, clearly 
showing metallic behaviori^iiSiii. 

In addition to that, in an extensive numerical analysis 
of the two-dimensional Anderson model with dimerised 
disorder, we have reported the existence of several dy- 
namical regimesi^ii^. 

As far as we are concerned there is no a theoretical ex- 
planation that describes adequately the metal-insulator 
transition in two-dimensional systems, as well as the dra- 
matic increase of the spin susceptibility in its vicinity. 
Comprehensive studies of the state of the art of this in- 
triguing problem are presented by Abrahams et al^ and 
Kravchcnko and Sarachik---. 

In one-dimcnsional systems with random disorder of 
any intensity, all states are exponentially localized as 
shown in the pioneering work by Anderson-^ in dealing 
with diagonal disorder, i.e., a model where the on-site 
energies are randomly distributed. On the other hand, 
when some correlation is included in the model, and with- 
out considering interaction between particles, this picture 
is substantially modified, given place to the appearance 
of extended states and, consequently, carriers are able 
to propagate through the system. Several structures 
with correlated disorder show vanishing of localization, 
when one considers nearest neighbors hopping. Among 
them one can quote the random-dimer model that can 
explain the high conductivity of polymersi2ii2iS2i2i. An- 
other example that show correlated disorder, responsi- 
ble for particle diffusion in ID is provided by the struc- 
tures where the on-site energies follow the Fibonacci and 
Thue-Morse sequence32242iS^i2^i2&. One should mention 
also the Harper model of a quasicrystal which presents a 
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mobility edge when the strength of the potential equals 
the half-bandwidth. Starting with a well localized par- 
ticle in the lattice, as long as the Harper potential is 
less than half-bandwidth, we encounter ballistic propaga- 
tion 27.28,29,30^ rp-^_^^ purpose of this work is to analyze the 
role the electron-electron interaction plays on propaga- 
tion in some ID nonperiodic structures. First we present 
in Sec. II the model assumed for the interaction between 
two particles in a band, namely the Hubbard Hamilto- 
nian. The two electrons are assumed in the singlet state 
in order to detect the effect of the Hubbard term, since 
it acts on opposite spins. We show the density of states 
for different strengths of the interaction. The dynamics 
tools introduced in order to characterize the dynamical 
behavior are presented in Sec. HI, namely, the time evo- 
lution of the mean square displacement (MSD) and the 
centroids associated with each of the particles, as well as 
the construction of 3D graphs of the wave packets evo- 
lution. In sec. IV we discuss the interplay between the 
strength of the interaction U with the initial wave packet 
extension. In Sec. V we present the results concerning 
to a ID crystal with impurities. In Sec. VI we present 
the conclusion to which we arrived in this work. 

I 



were = £r + Fedr, Sr being the on-site energy, F is 
the intensity of a dc electric field, c+^(cr.,s) is the Fermi 
creation (destruction) operator for an electron of spin 
s at site r, V is the hopping term and nrT(^rx) is the 
number operator for spin up (spin down) at site r. As it 
was said above, in order to analyze the role the U term 
plays on the dynamics, we treat the case of electrons with 
opposite spins, the singlet. 

A. Energy spectrum for a crystal 

To obtain the energy spectrum for the singlet in 
an impurity-free ID lattice, we solve the stationary 
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In this equation the first index refers to a particle with 
spin up and the second for spin down. The Wannier 
amplitudes 5nT,mido not depend on time. For simplicity 
we omit the label E in the Wannier amplitudes. 



II. THE HUBBARD HAMILTONIAN FOR TWO 
ELECTRONS INTERACTING IN A BAND 



With the aim to study the influence of the Coulomb 
interaction between carriers in a 11? lattice, we treat the 
electron-electron interaction along the Hubbard model^i 
in a crystalline system with impurities. As the ex- 
periments show, the electron-electron interaction should 
modify the behavior of the carriers as obtained in a one 
particle non-interacting scheme. In our work we intro- 
duce two interacting particles in an otherwise empty 
band. The study of such a problem have deserved a 
number of interesting works due to the relevant question 
of what is the role of the interaction between electrons 
on propagating properties in low dimensional disordered 
systema^^i^^i^'*i^^i^^i^'''i^^ . We assume that by solving the 
present problem one can get a better understanding of 
the role the interaction plays in real systems. 

We consider a one-dimensional lattice of N sites with 
lattice parameter d, for which the Hubbard Hamiltonian 
is: 



(1) 
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Schrodinger equation in the Wannier representation 
where we expand the eigenfunction in terms of the kets 
\ns, msi > that represent the state with one electron of 
spin s at site n and the other with spin si at site to: 



= g{ns,'ms'; E) \ ns,ms' > (2) 

ns,ms' 



In the Wannier representation we obtain the following 
set of equations corresponding to energy E: 



(3) 



For the crystalline case and without the presence of an 
electric field, all can be taken as zero. Introducing the 
center of mass and relative coordinates in units of the 



H 



+ C.C.) + U 2, flr^hri 



^ {dn+l.m + .9n,m+l + .9n-l,m + 3n,m-l) + {jn + Im + U 5n,m) gn,m — E gn,m 
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lattice parameter: 



R = (n + m)/2; 



r = n — m 



(4) 



and following Hubbard we expand the Wannier ampli- 
tudes: 



the electron-electron interaction in a crystalline structure 
in ID have produced a band with a structure similar to 
a 2D band for one electron, plus an extra band provided 
U is great enough. 



III. THE DYNAMICS 



(5) 



to obtain the following equation for the ^{K, k) : 



We follow the time evolution of wave packets in a ID 
lattice, by assuming the following expansion for the wave 
function in the Wannier representation 



E - 2y cos(f + k)- 2ycos(f - k) 
which in turn becomes: 

U ^ 1 



(6) 



(7) 



For every K we obtain — 1 roots of this equa- 
tion inside the infinities of the right hand side for which 
E— eK^f.+ei<_i. — 4Fcos-yCosfc. At the same time, 
and provided U is large enough, there appears an extra 
root outside the band. In this case, by varying K we 
get an excited second band. By looking at the eigen- 
value equation Q we note its close resemblance with the 
corresponding Schrodinger equation for one particle in a 
tight-binding 2D Hamiltonian, except for the U Hubbard 
terni'^^. The obtained density of states for the crystalline 
case for U ~ 5V shown in Fig. ^a, present two bands. 
The lower band is very similar to the case of one elec- 
tron in 2D for the reasons stated above and, the excited 
band has the typical structure of tight-binding ID bands. 
Moreover, by looking at the equation {Tj) one notes that 
E = U is a root of it, which happens for K — ±7r. This 
way if U is less than 4V we obtain resonant states inside 
the band, instead of an excited band separated from the 
2D band for a U—AV gap, that appears when U is greater 
than AV, since the bottom of the excited band occurs at 
E = U. See Fig. ^b. We conclude that the inclusion of 
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FIG. 1: The density of states in a crystal for: a) C/ = 5, and 
b) (7 = 1. The energy is in units of the hopping parameter V. 



where the Wannier amplitudes now depend on time. The 
time evolution of the wave function 



^-L^ - iJ I * > 

ot 



becomes in the Wannier representation: 



(9) 



J 



ih 



dfn 



dt 



^ (/n+l,m fn.rn-\-l /■n — l,m H~ fn.r> 

I 



l) + {ln+ 7m + U5n,-rn) fn 



(10) 



We assume the initial condition with an electron 
around site p with spin up and the other around site 
q with spin down. In order to do this, we shall consider 
the initial condition described by gaussian packets of dif- 



ferent standard deviations centered at sites p and q: 



I 4'(t = 0) > = |pT,q^ > 



(11) 
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where 



< n I p| (0) > = C exp 



< m I qi(0) >= Cexp 



(n - pf 



2a2 



(to — qY 
20-2 



(12) 



(13) 



and C is a normalization constant. 

As said before, all 7^ can be taken as zero for the crys- 
talline case without an electric field. We shall consider 
gaussians with standard deviations a. The case ct = 
corresponds to a particle localized in a single site. In the 
crystalline case where all are the same and not taking 
into account the Hubbard electron-electron interaction, 
both particles should propagate ballistically in the reg- 
ular lattice^. But now, in the presence of the on-site 
interaction C7, a different behavior should take place. To 
follow the evolution of the particles injected around sites 
p and q with spins up and down respectively, we perform 
the following. We define the mean-square displacement 
of each electron as indicated below. 



IV. INTERPLAY BETWEEN U AND THE 
INITIAL WAVE PACKET EXTENSION FOR A 
CRYSTALLINE CASE 

In order to study the combined effect that the interac- 
tion and the extension of the initial wave packet have on 
the propagating properties, we shall analyze the behavior 
of the MSD's as function of time. We consider the case 
in which both particles are injected at t ~ in the same 
position in the lattice. By taking four different initial 
waves associated with cr = 0, 1, 2 and 3 we note ballistic 
propagation in all cases, i.e., MSD = Dt^ + C. In Fig. [21 
we have plotted the ballistic coefficient D as function of 
the interacting Hubbard term U. 

We treat first the case of non-interacting electrons, i.e., 
U = We note that the MSD's values decrease when 
the extension of the initial packets increase. The broader 
the packet the more inertia it has, it becomes more "mas- 
sive" . 

We consider now the interacting case. First, we con- 
centrate on the CT = case and vary the strength of the 
Hubbard potential. The U term acts as an on-site en- 
ergy so that by increasing it, we introduce a barrier that 
inhibits hopping, the more so, the bigger U is. We still 
have propagation but with smaller MSD's values as we 
increase U. The ballistic coeflacient D is a monotonically 
decreasing function of U. 



and 



MSD,^ = J2(m-qf I 



(15) 



At the same time we can follow the centroid of the wave 
packet associated with the movement of each electron as 
follows: 



< Ax{t) = ^ (ri - p) 5] I fn,mit) I' (16) 



and 



< Ax{t) >ql = X! ("* ~ H I /",'"(*) 



(17) 



which will give us the amount of the displacement from 
the initial positions for both electrons. 

At the same time, by taking the absolute value of the 
difference between the respective centroid positions, we 
can determine the time evolution of the mean distance 
of the two electrons, or the size of the propagating pair: 
d{t) = \ <x{t) >pt - <xit) >gi\. 

In what follows we shall define the energies in units of 
the hopping parameter and the time in dimensionless 
units r = Vt/U. 




U (in units of V) 

FIG. 2; The MSD coefficient D as function of the interaction 
U for four cases of initial wave packets: cr = 0, 1, 2 and 3. As 
said in the text MSD/, = /, D/, + C 

An interesting situation occurs when the initial wave 
packet has an extension, i.e., cr 7^ 0. At this point we 
should state that the hopping parameter was assumed 
to be negative. By looking at the dynamical equation 
(10) we realize that what is relevant in analyzing prop- 
agation is the difference U - V, which can be considered 
like an effective interaction. This way, when both mag- 
nitudes have the same sign the interaction is reduced, 
while when they have different sign the effective inter- 
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action is increased. Consider first the case U > 0. We 
obtain in all cases studied {a — 1, 2 and 3) the D coef- 
ficient increases rapidly with U until reaching a plateau 
while D{(Ti) > D{a2) > D{a-i), i.e., the more extended 
the initial wave packet, the more massive it is. This also 
happens for the non interacting case J7 = 0, as shown 
above. 

Now for [/ < a different behavior is observed, where 
in this case the effective interaction is reduced. First 
we take a — I and see that for U = —1 smaller values 
than for U — are obtained, while for U= -2 we get even 
smaller values than for [/ = — 1. From then on the MSB's 
increase monotonically with \ U \, such that from | ?7 | = 
5 on, the values are greater than the corresponding to 
U = 0. The minimum of the coefficient D occurs for U = 
— 1.5. Taking a more extended initial wave packet, a — 2, 
the same trend is observed, but now the minimum is at 
U = -0.75. Aheady for = -2 the MSB is larger than 
for U — 0, and from then on they increase monotonically 
with I [/ I , as it happens for a ~ I. As for (t = 3 it shows 
the same behavior with the minimum closer to U = 0. 

We will explain now the initial decreasing trend for 
small \ U \ , when the initial packet has extension different 
from zero and hopping and interaction parameters have 
same sign. We have two competing effects to consider, on 
one side the extension of the packet implies more inertia 
(smaller MSB values) as said above, on the other hand, 
the presence of the U repulsive term acts on the different 
components of the wave trying to push the packet further 
away. For small \ U \ , the first effect overcome so one gets 
smaller MSB values but, increasing \ U \, the repulsion 
dominates so greater values are obtained. Finally, we 
note a cross over region that happens for 4 < | U \ < 5, 
where for the three cases considered: a = 1, 2 and 3, the 
corresponding MSB values are greater than the ones for 
cr = 0, for 1 [/ I > 5. 



V. TWO ELECTRONS IN A CRYSTAL WITH 
IMPURITIES 

We shall consider injected particles at time zero around 
sites p and q, where these sites as well as their nearest 
neighbors are taken as impurity levels in an otherwise 
regular lattice. In Fig.|21 we show the initial wave packets 
for the configurations ct = and 1. We illustrate the cases 
in which the packets are centered at positions ±2, and 
the impurity levels at e = —2, see Fig. |31 Note that for 
the case a — 1 the initial wave packets overlap due to 
the vicinity of the sites p and q. Obviously for cr = no 
overlap is obtained. 
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FIG. 4: Show for e = -21/, = 4 and cr = 0, in the up- 
per part, the centroids and the MSD as functions of time 
and in the lower part the time evolution of the up and down 
spin wave packets. Note that a considerable amount of both 
waves get trapped around the starting positions, while the 
rest propagate rapidly given rise to the observed superballis- 
tic behavior. 
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FIG. 3: The initial wave packets for the two configurations 
assumed: well localized (cr = 0) and extended (cr = 1). 



First we assumed the impurity energy levels at e = 
—2 and varied cr as well as the interacting factor U. We 
shall discuss the results for the initially well localized 
wave packets (cr = 0) and moderate interaction, J7 = 4, 
see fig. 0] The centroids clearly start to depart from 
each other as soon as the up and down packets overlap 
due to the interaction. The MSB shows a superballistic 
behavior: Ci" with a time exponent a = 2.05. One notes 
by looking at the wave packets evolution that rather big 
components of the two waves get trapped around the 
starting positions due to the well depth. At the same 
time, the tails of the initial packets are able to propagate, 
being responsible for the superballistic behavior. 

Now by taken cr = 1 we obtain a completely different 
picture, we observe an oscillatory movement of the cen- 
troids, see Fig.El The wave corresponding to each of the 
particles get trapped in a very small region and one no- 
tices that they strongly overlap in such a way that when 
one is moving to the right the other perform a displace- 
ment to the left and vice versa. The up and down spin 
wave packets get entangled with each other. The trap- 
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FIG. 5: The same as in Fig.Hbut for e = -2V, [/ = 4 and cr = 
1. In this case the more extended initial wave packets overlap 
given rise as a result to an entanglement of the two particles 
that perform oscillatory movements as time goes. Note that 
the MSD values are much smaller than in the previous case 
cr = 0. 



FIG. 7: The same as in Fig. H but for e = -2V, i7 = 8 and 
(7 = 1. Note the analogous behavior as the shown in the case 
for [7 = 4. 



the centroids is larger in the present case, see Fig. [7| 



ping of the wave packets can be understood by noticing 
that the sites covered by the gaussians are degenerate 
with on-site energies very different from the rest of the 
lattice, thus inhibiting hopping to distant sites. The be- 
havior of the MSB's proportional to is due to the es- 
caping of the tails of the initial wave packets that can hop 
to the near degenerate sites. This case serve to indicate 
that the propagating properties can not be inferred only 
from the MSD. In this connection we resort to 3D graphs 
that show the way the packets evolve in time. 






FIG. 8: The same as in Fig. H but for e 
a = 0. 




FIG. 6: The same as in Fig. H but for e = -2V, (7 = 8 and 

(7 = 



Increasing the interaction factor to [/ = 8 and taken 
cr = we obtain a similar behavior than the one described 
above for U — 4, see Fig. El For the same U and cr = 1 
we observe the same entanglement that occurs for J7 = 4, 
with the difference that the period of the oscillation of 



We treat now the case in which the impurity levels are 
deeper, i.e., at e = —4. For the initially well localized 
packets, cr = 0, and the interaction U = 4, we note that 
the centroids tend to get near to each other until they 
fill the repulsion, but they do not depart very much and 
later on they repeat the tendency to get closer. As for the 
MSD's the values are much smaller than the correspond- 
ing for e — ~2. By looking at the evolution of the wave 
packets in the 3D graphs it is evident that the packets get 
trapped due to the deepness of the wells, see Fig.|Sl Con- 
sidering now more extended initial wave packets, cr = 1, 
we note that the centroids depart from each other, the 
MSD's show ballistic behavior due to the escaping tales 
of the initial wave packets. As for the time evolution of 
the up and down wave functions we observe that they get 
trapped around the starting positions. A small amount 
of the amplitude of each of the up and down waves show 
a tendency to overlap with the other, see Fig. |51 
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FIG. 9: The same as in Fig. H but for e = -W , )7 = 4 and 
a = 1. 



Now for a strong interaction, [7 = 8 and taken ct = 0, 
we observe that the centroids tend to approximate while 
the MSB's values are much smaller than the correspond- 
ing to £ = —2. Due to the deeper impurities levels, the 
packets get trapped in a very small region of the lattice, 
see Fig. ^| When we take more extended initial packets, 
i.e., (7 = 1, the centroids tend to depart from each other, 
the MSB's are larger than the corresponding to the for- 
mer case, (7 = 0, and the wave packets perform the same 
kind of oscillations, see Fig. [TTl 

The time limit taken in our calculation was lO^^^sec, 
longer than any reasonable collision time, this implies 
that in order to avoid undesirable boundary effects, we 
have included lattices up to 1000 sites. We used the 
Runge-Kutta method to integrate the equations of mo- 
tion. 



VI. CONCLUSIONS 
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FIG. 10: The same as in Fig.|I|but fore = ~AV , t/ = 8 and 
a = 0. 
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FIG. 11: The same as in Fig.|I|but fore = ~AV , i7 = 8 and 
a = 1. 



In this work we analyzed the role the electron-electron 
interaction plays in one-dimensional disordered systems, 
in particular, a crystalline lattice with impurities. Two 
electrons with opposite spins were introduced in a band 
where their interaction was assumed along the Hub- 
bard model. We have obtained the energy spectrum of 
a crystal considering different values of the interaction 
strength. The obtained density of states show that for 
sufhciently large C/, two bands are present separated by 
a gap of intensity U — AV . The lower band has the fea- 
ture of the one corresponding to a single particle in a 2B 
crystal, while the excited band has the structure of a IB 
tight-binding crystal. 

As for the dynamical properties, we analyzed the in- 
terplay between the initial wave packet extensions and, 
for both, positive and negative values of the interaction 
U. 

In the non-interacting case, U = 0, we show that by 
increasing the extension of the packets, one reduces the 
MSB values, the more so, the more extended the initial 
wave packets, they become more "massive". 

In the case of a positive U, the obtained MSB increases 
rapidly from the non interacting case [/ = 0. For neg- 
ative [/, a different initial behavior is observed. As U 
is near the origin the MSB decrease until a minimum 
value is achieved, and from then on it increases reaching 
a plateau. For both the positive and negative values of U 
there exists a cross over region as explained in the text. 

Next we analyze the propagation that takes place when 
we introduce the two interacting particles at impurity 
sites in an otherwise periodic crystal. We have encoun- 
tered situations for which the MSB as function of time 
presents superballistic behavior, although when looking 
at the packet we realize that most of it get trapped in 
a small region of the lattice. That is why the propagat- 
ing properties can not be inferred from the MSB alone. 
When we consider moderate impurity levels, in our case 
£ = —2V, the response of the system with regard to the 
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initial conditions arc very different. In fact, when we 
start with very well localized particles, ct = 0, the cen- 
troids depart from each other for every intensity of the 
interaction, while for the less localized initial packets, 
a = 1, the two particles get entangled, trapped in the 
impurity region, we observe the centroids performing a 



periodic movement. When considering deeper impurity 
levels such as e = —4 no matter what the extension of 
the initial waves packets are, the result indicate that the 
particles remain localized around the starting positions, 
this is true for every interaction strength. 
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